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ABSTRACT 


The  equations  of  motion  of  a  body  under  the 
influence  of  the  earth' s  gravity  and  atmospheric  drag 
are  obtained  in  radar  polar  coordinates.  The  choice 
of  this  coordinate  system  may  have  important  advantages 
in  various  tracking  and  filtering  algorithms.  A  rotating 
ellipsoidal  earth  with  gravitational  terms  up  through  the 
second  (easily  extendable)  harmonic  is  assumed.  The 
effects  of  a  non -spherical  earth  along  with  the  corres¬ 
ponding  effects  of  the  higner  gravitational  harmonics 
can  be  important  in  the  estimation  of  ballistic  coefficients 
at  verv  low  deceleration  levels. 
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GLOSSARY 


r  vector  from  earth  center  to  target 

R  earth  radius  vector  to  "foot"  of  radar  site 

S  vector  from  radar  to  target 

S,  £,a  radar  range,  elevation,  and  azimuth,  respectively 

s,£,a  unit  orthogonal  vectors  in  the  directions  of  increasing  S  , 
e,  and  a  ,  respectively 

X,  Y,  2  rectangular  coordinates  with  origin  at  radar  and  X  East, 

Y  North,  and  Z  vertically  up 

1 ,  j\  k  unit  vectors  along  X,  Y,  2,  respectively 

c  unit  vector  along  earth  axis  £  ,  pointing  toward  North  Pole 

r  unit  vector  from  earth  center  to  target 

<ju  earth  rotation  rate 

A 

X  earth  rotation  vector  (uT=  u>  §) 

u  geodetic  latitude  of  radar 

pc  geocentric  latitude  of  radar 

cr  geocentric  latitude  of  target 

R  earth  radius  at  latitude  cc 

<r 

a  equatorial  radius  of  earth  (6,  378.  160  km) 

H  height  of  radar  above  ellipsoidal  earth  (i.  e.  ,  above  MSL) 

I#  moment  of  inertia  of  earth  about  its  axis  £ 

I  moment  of  inertia  of  earth  about  an  axis  in  equatorial  plane 

GM  product  of  gravitational  constant  and  earth  mass 

(3.  986032  x  105  km3/*ec2) 

A  gravitational  acceleration  vector  at  the  taiget  position 
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drag  velocity  vector  (i.e.  ,  relative  to  air) 
inertial  acceleration  vector 
drag  acceleration  vector 

height  oC  target  above  ellipsoidal  earth  (i.e.  ,  above  MSL) 
atmospheric  mass  density  at  target  position  (p  =  ?(h)  ) 
acceleration  of  gravity  at  target  position  (g  =  |A  |) 
mass  of  target 

ballistic  coefficient  =  m/CjyA  (  i.e.,  mass-to-drag  ratio) 

eccentricity  of  reference  ellipsoidai  earth  (e^  =  .0066945) 

2  -3 

dimensionless  constant  =  3  (Ip  -  Ig)/(2Ma  )  =  1.624  x  10 
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Trajectory  Equations  of  Motion  in  Radar  Poiar  Coordinates 
INTRODUCTION 

Trajectory  equations  of  motion  in  radar  rectangular  coordinates  (the 
X,  Y,  Z  system  below)  are  readily  available  in  the  literature.  They  are  not 
to  the  author' s  knowledge  available  in  radar  polar  coordinates  (Range  S, 
Elevation  e  ,  and  Azimuth  a  }.  Perhaps  because  of  this  (and  possibly  be¬ 
cause  of  the  fear  that  such  equations  would  be  too  complex  anyway),  it  has 
been  common  practice  in  problems  of  radar  tracking  to  do  the  estimation  and 
prediction  in  radar  rectangular  coordinates  even  though  the  data  are  gathered 
and  the  radar  is  directed  in  radar  polar  coordinates.  This  procedure  involves 
wastefuL  coordinate  transformations  in  at  least  one  direction  which  can  be 

avoided  completely  if  radar  polar  coordinates  are  used  throughout.  These 

♦ 

coordinate  transformations  are  particularly  complex  if  range  rate  (S)  and 
range  acceleration  (S)  are  used.  On  the  other  hand,  the  equations  of  motion 
in  radar  polar  coordinates,  derived  in  this  report,  turn  out  to  be  only  moder¬ 
ately  more  complex  than  in  radar  rectangular  coordinates.  Quantitative  com¬ 
parison  of  the  computational  economy  of  using  one  coordinate  system  or 
another  will  depend  on  the  details  of  the  specific  application.  For  example, 
if  the  data  rate  is  high  and  large  step  sizes  can  be  tolerated  in  the  integration 
of  the  equations  of  motion,  it  may  be  more  economical  to  use  radar  polar 
coordinates  throughout,  while  in  the  reverse  situation,  the  use  of  radar 
rectangular  coordinates  may  be  advantageous.  Additionally,  the  following 
important  advantages  of  radar  poiar  coordinates  in  connection  with  track 
initiation  in  Kalman  filtering  have  been  brought  out  in  discussions  with 

H.  M.  Joisc-s  and  R.  P.  Wishner  of  Lincoln  Laboratory: 

I.  The  covariance  matrix  of  the  estimation  parameters  in  the 
Kalman  filter  is  more  nearly  diagonal  in  radar  polar  co¬ 
ordinates  and  thereby  less  susceptible  to  the  effects  of  ill- 
conditioning  with  short  lengths  of  track 

2.  The  usual  linearity  assumptions  of  Kalman  filters  are  more 
nearly  correct  in  radar  polar  coordinates.  This  can  result 
in  better  track  initiation  since  initial  estimates  can  be  less 
accurate  and  still  yield  a  nearly  optimum  correction  with  later  data. 
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DERIVATION 

Assuming  that  only  gravity  and  drag  forces  act  on  the  body  of  mass,  m  , 
we  have  for  the  equation  of  motion 

F  =  F  ,  +  F 
d  g 


or 


mA  “  mA  ,  +  mA 
o  d  g 


The  drag  force  is  taken  to  be 


(1) 


Fd  =  m\i  =  •  <2> 

,  # 
where  (3  =•  m/CjyA  is  taken  as  the  definition  of  ballistic  coefficient  and  o 

is  the  atrnr  f’~ric  mass  density  at  the  body  position.  Using  Eq.  (2)  in 

Eq.  (1)  and  div.  iing  by  m  gives  for  the  equation  of  motion 


Ao  =  *'#V+Ag  •  O) 

V  is  the  target  velocity  relative  to  the  moving  radar  site  ior  equivalently, 
relative  to  the  surrounding  air,  assuming  the  atmosphere  to  rotate  rigidly  with 
the  earth  and  no  local  winds)  and  can  be  immediately  expressed  in  the  §,  e,  a 
system  as 

V=Vs+Ve+Vd,  (4) 

sen 

where 


V  =  Se 

£ 


V  ~  Sa  cos  € 


*  This  definition  avoids  the  needless  confusion  which  the  definition 

{3  =  mg/Cj-jA  introduces  as  to  whether  tc  evaluate  g  at  sea  level  or  at  the 

target. 
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A  is  the  gravitational  acceleration  vector  and,  following  McCuskey, 
§> 

can  be  expressed  as 


.  A  /> 

A  = g  r  + g_^ 
g  &r  B§' 


where 


gr  =  _  t 1  +  ‘T  (7) 2  ( 1  -  5  sin Z(p)  3 

r 

2G  M  r  ,  a ,  2 

g*  - - ”  J  (r}  n  V  * 

b  r 


Its  magnitude  is  given  by 

g  =  |Ag{  =  [Ag  •  Ag ] ^  =  [gr2  +  g^2  +  2gr  g?  »in  (p^Z  . 

In  order  to  express  Eq.  (5)  in  the  s,  e,  a  system,  it  is  necessary  to  express 
r  and  Z  in  that  system.  This  is  outlined  in  the  next  paragraph. 

The  geocentric  target  position  vector  r~  is  given  by  (see  Fig.  1  and 
glossary) 


Since 


r=R+Hk+S 


A  *■  rt 


S  =  Xi  +  Yj  +  Zk 


Nsg,- 

f-*£?r  - 

1  W’- 

j  ;T 


we  nave 


where 


R  -  -  R  sin  (u  -  u  )  j  +  R  cos  (u  -  u  )  k  , 
c  « 

r=  Lx]  i  +  [Y  -  RQ2!  j  +  [Z  +  RQa]  k 

Qj  =  cos  (n  -  pc)  +  H/R 

Qv  -  sin  (p  -  u  ) 

tf  v. 


*•  $.  W.  McCuskey,  ’’Introduction  to  Celestial  Mechanics, ”  pp.  163-166  (1962). 
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I  -  e  cos  y 


tan"1  [(1  -  e  )  tan  u] 


AAA 


Using  the  transformation  (12),  we  obtain  r  in  the  s,  £  ,  a  system  as 


-r  ^  A  ^ 

r  -  r  s  +  re  +  r  a 


where 


r  =  S  +Q,R  sin  e  -  Q,R  cos  e  cos  a 
si  c 

r  =  Q,R  cos  e  +  Q,R  sin  e  cos  a 

cl  *■ 


r  =  Q-,R  sin  a 
a  2 


Its  magnitude  is  given  by  j 

r  =  j7|=  Rl(S/R)Z  +  2(S/R)(Q1sine -Q2cosecosa)+Q1  +QZ  } 


.  (7) 


The  unit  earth  rotation  vector  %  can  be  expressed  as 


A  A  .  * 

§  =  cos  pj  +  Sin  ilk 


Again,  utilizing  the  transformation  (12),  this  can  be  expressed  in  the 
system  as 


AAA 

s,  e,  a 


where 


A  a  _  A  a 

E=Es+Ee+£a 

"CL 


c  =  cos  ii  cos  c  cos  Q  +  sin  n  sin  c 
■"S 

^  =  -  cos  (i  sin  e  cos  a  +  sin  ^  cos  e 

p  =  -  cos  a  sin  a 
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The  required  quantity  sin  <p  is  then  simply  obtained  as 

— ♦  A 

sirup  =  =—  {S  cos  e  cos  a  cos  p  +(H  +  S  sin  e)  sin  p  +  R  sin 

Using  Eqs.  (6)  and  (8)  in  (5)  yields  A  in  the  s,  e,  3  system  as 

o 


—  e  r 
Br  s 


,grr 


g  r 

°r  a 


A  +  gcE  )s  +  (- t-2-  +  g*  )  e  +(—=——  +  g„F  1 

p  '  r  bC's'  v  r  &?-£'  '  r  ®c'Q 


(9) 


(10) 


Equations  (7)  and  (9)  are  used  to  evaluate  and  .  All  that  remains  is  to 

express  Aq  in  the  s,  e,  a  system  which  follows. 

A  ,  the  earth-centered  inertial  acceleration  of  an  object  observed  from 

$ 

a  site  on  a  uniformly  rotating  earth,  can  be  expressed  as 


A  -(V  s +V  £  +  V  a +V  il +V  £ +V  &)+2ujxV +uTx(uTxr) 
o  seassa 

where  the  term  in  parentheses  (V)  is  the  apparent  (i.  e.  ,  relative  to  the 
moving  radar  site)  acceleration  of  the  body. 

The  cartesian  and  polar  unit  vectors  are  related  according  to  the 
transformation 


r  A  H 

S 

- 

I 

A 

"  i 

r*A 

s 

A 

€ 

=  M 

A 

j 

or 

A 

j 

=  m"1 

s 

A 

a 

w  J 

-  « 

A 

k 

A 

a 

»  m 

(ID 


(12) 


where 


cos  s  sin  a 
M  =  -  sin  G  sin  a 

cos  a 

.-1  -„T 


cos  e  cos  a  sin  €  * 

-  sin  e  cos  a  cos  g 

-  sin  a  0 


(13) 


and  M  =  M  (i.  e.  ,  transpose  of  M)  .  The  unit  vector  s  s 
lated  to  their  time  rates  as 


A  A  A 

l,  a  are  re- 


*  L..  Page,  "Introduction  to  Theoretical  Physics,"  3*u  Ed.  ,  4th  Printing, 
p.  121  (1959). 
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+  a  cos  e' 


-  a  sin  e 


a  cos  e  a  sin  e  0 


The  relation  (14)  is  easily  verified  by  differentiating  Eq.  (12)  with 

^  A  A  A 

“i  -  j  =  k  =  0  and  solving  for  s  ,  e  a  .  It  can  also  be  verified  directly  by  re¬ 
ferring  to  Fig.  1  and  using  relations  of  the  type 


a  ds  ds  Sc  ,  dS  da  3s 


dt  de  dt  +  da  dt  *  "§*£ 


=  e  ,  and  -Jr*  =  cos  e  a 


Forming  the  dot  product  of  Eq.  (11)  with  s,  a,  successively,  yields 

A  "  s  =  [V  +(V  e+V  4).  S]+2(aTxV)-  s+[uTx(iTxrl  ]*  s 
o  s  e  a  ' 

Aq*  c  =  [Ve  +  (Vg§  +V^)  •  e  ]+2(uuxV)  *  e  +[iox(uuxr)  ]♦  e 

A  •  a  =  [V  +(V  s  +  V  €)  •  a  ]+2(U)xV)  •  a  +  [u)x(usc  r)  ]  •  a  .  (15) 

O  s  c 

—  A 

Using  a- =  U)5  and  Eqs.  (4),  (6),  (8),  (9),  and  ( 14)in  ( 15),  and  reducing  to 
scalars  yields 

A  •  s  =  [S -S(fe2+a2cos2G)]  +  2u)[§  V  -  5 V  ]  +  [ Z  r  sin <p  -  r  ] 

A  •  e  =  [Se +2Se +Sa2sinecos  ej  +  2u)[5  V  -5  V  ]+u)2[e  rsixup-r  ] 

A  >  a  =  [Sacos  e  +  2(Sa  cos  e  -  Sae  sine)  ]+ 2u)[§  V  -  E  V  ]+u)2[§  r  sin</>  -  r  ] 
o  sees  *  »* 

Noting  that  ^  ^ 

Aq  =  (Aq  •  s)  s  +  (Aq  *  e)  e  +  (Aq  •  a)  a  (17) 

we  finally  have  Eq.  (3)  expressed  completely  in  the  s,  e,  a  system  by  using 
Eqs.  (4),  (10),  and  (17).  Expressing  Eq.  (3)  in  scalar  form,  we  obtain  the  set 
of  three  differential  equations  of  motion  in  radar  polar  coordinates.  The 
result  is 


-H*  W'lWtlWlllw. 


I 


S=S(e2  +  a2cos2e)  -  2u<§eVQ-  ?QV£ 


V  =^-{2Se+Sa2sinecose+2u){§QVs-?sVQ)+u)‘'{5ersin</?-re)  +  (-^)Ve 


5„V„-  §nVc)  -w2(5ar  sin?  -  rg)  -(^)Vs+£r^  +g£F5 


g  r 

—  -  g  §  } 


S  cos  e 


where 


{ 2(Sa  cos  e  -  See  sin  c)  +  2 u)(SsV£  -  ievg)  +«M§Qr  sin?  -  rG)+(|j£)  VQ 

g  r 

•  11 


f  =  cos  u  cos  e  cos  a  +  sin  p  sin  e 
-s 

£  =  -  cos  (a  sin  e  cos  a  +  sin  p  cos  e 


=  -  cos  p  sin  a 

=  S,  V  =  Se,  V  =  S h  cos  e 
6  “ 


=  (v  2  +  v  2  +  v  2}V* 

s  e  a. 

=  S  +  Q  jR  sin  e  -  Q2R  cos  e  cos  a 
=  Q  jR  cos  e  +  Q^R  sin  e  cos  a 
=  Q^R  sin  a 
=  cos  (p  -  pc)  +  H/R 
=  sin  (p  -  pc) 


=  R  {(S/R)2+2(S/R)(Q1sir  e  -  Q^osecosa)  +  Q^+Q^}1^ 


ll  +  J  (f)2  (1-5  Sin2?)} 
r 

J(j)2sin? 

r 


sin? 


=  i.  {s  cos  e  cos  a  cos  p+(H+S  sin  e)  sinp+Rsinp  } 

1*  C 
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The  height  h  required  in  p  =  p(h)  is  obtained  as  the  distance  between  the 
target  and  the  point  of  intersection  of  the  reference  ellipsoid  and  the  line 
passing  through  the  target  and  normal  to  the  reference  ellipsoid.  An  ex¬ 
cellent  approximation  to  this  is 

1  I/* 


h  ^  r  -  R 


a2{l 


e2) 


1 


e  cos  cpj 


This  if-  the  distance  between  the  target  and  the  intersection  of  r  with  the 
reference  ellipsoid. 

These  equations  are  only  moderately  more  complex  than  the  equations 

of  motion  in  radar  rectangular  coordinates.  Only  the  first  two  terms  in 

each  equation  are  unique  to  the  polar  coordinate  system.  Coriolis  accelera- 

2 

tion  terms  (containing  w  ),  centrifugal  acceleration  terms  (containing  u,  )  , 
drag  terms  (containing  c/26  )  and  gravitational  terms  (containing  g^  and 
g_  )  appear  in  both  sets  of  equations. 

The  relation  (3),  or  equivalently  the  relations  (18),  can  be  used  to  ob¬ 
tain  various  expressions  which  yield  estimates  of  P  from  radar  data.  Also, 
it  is  easy  to  specialize  to  the  case  of  a  spherical  earth  with  l/r^  gravity  by 
simply  setting  J  =  e  =  0  and  u  =  Pc  .  It  is  than  useful  to  observe  what  kind 
of  errors  can  result  in  p  due  to  earth  model  choice.  Computer  runs  based 
on  actual  radar  data  have  indicated  the  following.  The  ballistic  coefficient  , 

P  ,  estimated  by  using  the  simple  model  (J  =  e  =  0),  can  differ  from  that 

estimated  by  using  the  more  accurate  model  (J,  e  as  reported  here)  by  more 

_3 

than  a  factor  of  two  when  the  drag  acceleration  A  ,  is  g  x  10  .  The  radar 

U 

data  used  were  at  the  mid-latitudes  where  these  differences  are  probably 
greatest.  A  good  conservative  rule  would  be  to  require  use  of  an  "improved" 
model  when  attempting  to  estimate  ballistic  coefficients  at  drag  levels  in  the 
vicinity  of  one -thousandth  of  a  g  or  less.  An  "improved"  model  means  one 
at  least  as  good  as  that  reported  here  (i.  e.  ,  ellipsoidal  shape  with  gravitational 
terms  through  the  second  harmonics). 
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